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We systematically derive the consistency relations associated to the non-linearly realized symmetries 
of theories with spontaneously broken conformal symmetry but with a linearly-realized de Sitter 
subalgebra. These identities relate (A^+l)-point correlation functions with a soft external Goldstone 
to A^-point functions. These relations have direct implications for the recently proposed conformal 
mechanism for generating density perturbations in the early universe. We study the observational 
consequences, in particular a novel one-loop contribution to the four-point function, relevant for 
the stochastic scale-dependent bias and CMB /i-distortion. 



Although non-linearly realized symmetries have long been known to be a powerful tool in many 
areas of physics, only recently has this technology been brought to bear on theories of the early uni- 
verse. For example, it has been realized that perturbations in single-field inflation can be described 
most generally by the effective field theory of spontaneously broken time diffeomorphisms [1, 2]. 
Single-field inflation can also be understood in terms of global symmetries as the spontaneous break- 
ing of the 5*0(4, 1) conformal symmetry of down to spatial translations and rotations [3, 4]. The 
corresponding Goldstone field is the curvature perturbation of uniform-density hyper surf aces. 
Moreover, the well-known consistency relations [5-9], which constrain the soft limit of correlation 
functions, arise as Ward identities for the non-linearly realized symmetries [10-12]. Additionally, 
symmetry considerations have proven to be a powerful tool in analyzing correlation functions of 
spectator fields in inflation: both gravitons [13] and scalar field perturbations [14-18] are con- 
strained to have conformally- invariant correlators at late times. 

Recently, it has been realized that the de Sitter symmetries — including dilation invariance, 
responsible for the scale invariance of perturbations — need not correspond to the isometrics of a 
physical metric, as in infiation, but instead can arise as the unbroken sub-algebra of spontaneously 
broken conformal symmetry where the dynamical metric is nearly flat [19-21]. In this conformal 
mechanism, some fields in the conformal field theory (CFT) acquire specific time-dependent expec- 
tation values, which breaks the 50(4,2) conformal symmetry on (approximate) Minkowski space 
down to its de Sitter subgroup: 
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As a result of couplings dictated by conformal invariance, spectator fields in the theory evolve in 
a fictitious de Sitter background, and consequently massless fields acquire a nearly scale invari- 
ant spectrum of perturbations.^ The de Sitter expansion is fictitious; the physical, Einstein-frame 
metric describes a universe which is very slowly contracting or expanding. Incidentally, such slow 
evolution drives the universe to be increasingly flat, homogeneous and isotropic [23], thereby ad- 
dressing the well-known problems of standard big bang cosmology. A robust prediction of the 
scenario is the absence of gravitational waves — because the universe is approximately static, ten- 
sor modes are not appreciably excited. As with any mechanism relying on multiple fields, a generic 
prediction is a significant level of non-Gaussianity in the squeezed/local limit. 

An example of this mechanism is Galilean Genesis [20]. This scenario is based on the con- 
formal galileons [24], a class of conformally-invariant, higher-derivative scalar field theories, which 
nevertheless only propagate one degree of freedom. Galilean Genesis capitalizes on the fact that 
galileons can violate the Null Energy Condition (NEG) [25] in a stable manner, to describe a 
universe which expands from an asymptotically flat initial state. One drawback of the original 
scenario is that fluctuations around (slight deformations) of the NEG-violating background prop- 
agate superluminally. Superluminality can be avoided through explicit breaking of the conformal 
group, preserving only dilation [26]. Alternatively, the Dirac-Born-Infeld (DBI) generalization of 
the scenario [27, 28] preserves the full conformal symmetries, but, thanks to a different non-linear 
realization of the conformal algebra, features perturbations that propagate strictly sub-luminally 
around the NEG-violating background. Another incarnation of the conformal mechanism is the 
f7(l)-invariant model [19, 21, 29], which describes a complex scalar field rolling down a negative 
quartic potential. The resulting cosmology is a phase of slow contraction before a NEC-violating 
phase that leads to the big bang. 

In fact, the conformal mechanism is more general than these particular incarnations. All the im- 
portant physics — in particular scale invariance — is fixed by the symmetry breaking pattern (1.1), 
irrespective of the microphysical details. Applying the technology of the coset construction [30-32], 
the most general effective action describing the Goldstone field vr and other "matter" fields, includ- 
ing weight-zero fields, was derived in [33] systematically in powers of derivatives. By construction, 
the effective action linearly realizes the 5*0(4, 1) de Sitter algebra^ and non-linearly realizes the 
50(4, 2) algebra. Many key properties of the scenario follow immediately from this effective action, 
e.g., the scale invariance of massless fluctuations, and the fact that the time-dependent background 
is a dynamical attractor. 

In this paper we focus on the implications of the symmetry breaking pattern (1.1) on various 
correlation functions. As the models under consideration enjoy linearly-realized de Sitter symme- 
tries, we begin in Sec. 2 by reviewing some elementary properties of field theory on de Sitter space, 
focusing in particular on the late-time action of the de Sitter isometrics and their corresponding 
constraints on correlation functions. 

Our primary interest, however, lies in understanding how the non-linearly realized symmetries 
— time translation P^, boosts J^*, and the temporal component of the special conformal transfor- 
mations — act on correlation functions. These symmetries are discussed in Sec. 3. Analogously 
to single-field inflation, we will find that these result in consistency relations relating {N + l)-point 
correlation functions with a soft external vr to a (broken) symmetry transformation on the A^-point 

^It is well-known that weight-zero fields are forbidden by unitary bounds [22], but these assume a stable 
conformally invariant vacuum. The conformal mechanism does not assume such a vacuum, only a time- 
dependent symmetry-breaking background. The conformal vacuum can be unstable or not exist. 

^We will be somewhat loose about notation and refer to both Lie groups and their corresponding algebras 
in the same way, e.g., SO{A,2) for the conformal group/algebra. 
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Figure 1: Schematic representation of the consistency relations: {N + l)-point functions with a 
soft external vr leg are related to A^-point functions. 

functions without the soft vr. In Sec. 4, we derive the equal-time master consistency relation, which 
is the main result of the paper (see Fig. 1): 

\ a a / 

where the (/>'s are fields of arbitrary mass, the ' indicates a correlation function without the 
momentum-conserving delta function, and ^= indicates equality after performing an angular aver- 
age only over the terms quadratic in q. Since there is only one Goldstone for 5 broken symmetries, 
the above consistency relation constrains different powers of q in the soft limit g — )• 0. The 0{q^) 
part of the constraint results from , the 0{q^) part from j''*, and the 0{(f') part from K^. 
We derive (1.2) using the background- wave method [6], in which the long- wavelength vr mode is 
treated as a classical background for the short-wavelength modes; in a follow-up paper, we will 
show how (1.2) arises as a consequence of the Ward identities for the broken symmetries, analogous 
to the classic soft-pion theorems of chiral perturbation theory [34, 35]. We also check the validity 
of (1.2) in a myriad of cases using explicit correlation functions derived from the general effective 
action of [33] . 

The consistency relation (1.2) opens up the possibility of strong observational tests of the 
conformal mechanism, which is the focus of Sec. 5. In much the same way that observation of 
/nl in the squeezed configuration would rule out all single-field models of the early universe, 
observation of a violation of one of these consistency relations would rule out the production of 
density perturbations by the conformal mechanism. An important caveat in making contact with 
observations is the observability of the vr mode — although tt acquires a strongly red-tilted spectrum, 
its contribution to the adiabatic mode is actually strongly blue-tilted and hence is negligible on 
scales probed by cosmological observations. (The observed C, is dominated by the scale invariant 
contribution from conversion [36] of the weight-zero fields.) However, the consistency relation (1.2) 
does have observational implications for correlators in two ways. The first is for correlation functions 
with a soft internal line. In the limit where M external momenta sum up to a soft total momentum, 
the amplitude is dominated by soft-vr exchange, thanks to its strongly red-tilted spectrum, and 
factorizes into a product of (M + l)-point and {N — M + l)-point functions, each with a soft tt 
insertion: ^ 

^h, ■ ■ ■ hX^o = jr^i^-<i^k, ■ ■ ■ '^fcM)Uo(Mfc,,+i • • • hX^o ■ (1-3) 
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The consistency relation (1.2) then constrains each soft amplitude in the product. We show how 
this constrains the form of the four-point function for massless spectator fields x ™- the limit where 
the sum of two momenta is soft, reproducing the behavior explicitly calculated in [37-39]. We also 
find a novel contribution to the four-point function of spectator fields from a one-loop exchange of 
two vr's. Despite being a one-loop effect, it is enhanced compared to the tree-level amplitude by 
inverse powers of q so that it dominates in the soft limit. Its momentum dependence is of the tnl 
form 

TNL~^log|, (1.4) 

SO that we can expect signals in the stochasticity of the scale-dependent bias [40] and the power 
spectrum of CMB /i-distortion [41]. Although vr itself is not observable, there is still an observable 
signature of soft external vr fields; as was pointed out in [37-39], the presence of a long vr mode 
leads to statistical anisotropy in spectator correlation functions, in particular it leads to anisotropy 
of the power spectrum of x- 



l + ciV^^(3cos2 0-l)+C2p^cos2 01og^ ) , (1.5) 



where 9 is the angle between k and q, Hq is the present-day Hubble parameter, and ci^2 are constant 
coefficients which depend on the particular realization of the long- wavelength vr modes. We will 
show that this anisotropy follows solely from symmetry considerations, and is therefore also a 
generic prediction of the conformal mechanism. 

Since the scenario relies on de Sitter symmetries, it is natural to wonder whether the conformal 
symmetries can also be present in the inflationary context. The short answer is yes, but the 
inflationary realization is far less natural. Note first that the action we write down for vr cannot be 
the action for perturbations of the infiaton — they have linearly- realized 5*0(4, 1) symmetry, and 
thus do not fit into the effective field theory of inflation framework. If vr were the perturbations of 
the inflaton, the corresponding de Sitter space would be eternal, there is no clock telling inflation 
when to end. We are therefore led to consider a situation where vr is merely a spectator fleld^ 
in multi-field inflation — its presence therefore gives additional structure to the spectator sector 
and all of the same consistency relations will hold provided the entire spectator sector couples in a 
conformally invariant way to vr. However, this is a somewhat artiflcial situation. We can of course 
impose this additional structure, but it does not buy us anything new — whereas in the conformal 
mechanism this structure is a necessary and natural consequence of the mechanism. 

Various subtleties of our calculations are addressed in detail in Sec. 6, including complications 
introduced by the fact that vr is a fleld of negative 3d conformal weight in a de Sitter background, 
the ambiguity (or lack thereof) of off-shell versus on-shell correlators in the consistency relation, and 
the translation from the consistency relation (1.2) to the form expected from the Ward identities 
for spontaneously broken conformal symmetries. We summarize our results and discuss future 
spin-off directions in Sec. 7. Various appendices collect results peripheral, but important, to our 
main line of development: Appendix A derives the action of the late-time de Sitter isometries on 
correlation functions. Appendix B reviews the construction of actions for spontaneously broken 
conformal symmetry, and Appendix C computes the correlation functions needed to adequately 
verify the consistency relation. 



■^Of course, there is another alternative, namely for the field vr itself to acquire a profile and act as the 
inflaton. While possible, it is easy to check that the potential for vr required by conformal invariance only 
supports inflation over a limited region in field space, vr < Mpi, corresponding to approximately 1 e-fold of 
inflation. 
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2 Linearly realized SO(4,l) and 3d conformal transfor- 
mations 

The scenario relies on linearly realized S'0(4, 1) invariance, so we begin by reviewing some properties 
of scalar fields on de Sitter space. Throughout, we will work in the planar slicing of de Sitter space, 
where the line element takes the form^ 

= (-d*' + d^') ' (2-1) 

where t < is conformal time. This is identical to the situation in multi-field inflation, where 
spectator fields feel a de Sitter geometry and do not back-react appreciably. A key difference, 
as emphasized in the Introduction, is that de Sitter space is a fake geometry in the conformal 
mechanism — the actual, Einstein-frame metric is slowly evolving. Nevertheless, for the purpose of 
this discussion we can remain agnostic as to whether or not the background de Sitter corresponds 
to the actual metric. 

The de Sitter metric (2.1) corresponds to a maximally symmetric space-time and therefore 
enjoys 10 isometries. Six of these are the familiar translations and rotations of the flat spatial 
slices: 

x'^x' + a'- (2.2) 
x' J\x^ . (2.3) 

Additionally, de Sitter space is invariant under a dilation of both spatial and time coordinates 

x^" Ax^ . (2.4) 

Finally, it is invariant under the simultaneous transformation of space and time as 

t — >t-2t{h-x); (2.5) 
x' ^x' + h\-t^ + x^)-2x\b-x) , (2.6) 

where 6* is a real-valued 3-vector. 

Next, consider a free scalar field on the de Sitter background: 

S = /d^x^(^-i(50)2-^0^^ . (2.7) 

The de Sitter isometries act on (j) as follows: spatial rotations and translations (2.3) act in the usual 
way, 

5p^(j) = -di(j) ; 

^J^,(t' = {xidj - Xjdi)4> , (2.8) 
while the remaining four isometries (2.4) and (2.6) act as 

^D(t> = -{-tdt + X ■ d)(t) ; 

5kA = -(-2xitdt + 2xiX-d-{-t'^ + x'^)di)(P . (2.9) 



■^Here we use t as the conformal time coordinate on de Sitter space in order to emphasize the connection 
with models in which de Sitter arises as a fictitious background from broken conformal invariance. 
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We are interested in how these transformations act at late times (t — )• 0). In Fourier space, the 
equation of motion that follows from the above action in the coordinates (2.1) is 



2; /,2 



(t>k - -^k + \^k' + j^j (/.fc = , (2.10) 

with the well-known solution given by Hankel functions. In the long-wavelength {k — >• 0) limit, the 
time dependence of the mode functions simplifies to 



5fc 




, with A± = -±l ---f . (2.11) 



Assuming < 9H'^/4, the growing mode corresponds to A_ = A, and the time dependence of 
the field is (/> ~ as t — )• 0. In this limit, we can therefore replace tdt ^ A in the transformation 
rules (2.9) and neglect 0{t'^) terms to obtain 



6d(I) = i^A-x-dj 

= (2Axi-2xiX-d + x^di](l) . (2.12) 



These are recognized respectively as spatial dilations and special conformal transformations for a 
field of conformal weight A. Combined with the spatial rotations and translations (2.8), they form 
the conformal group on M^. Therefore, correlation functions of fields on de Sitter must be invariant 
under conformal transformations of Euclidean 3-space on the future boundary [13-17], which is of 
course the basis of the proposed dS/CFT correspondence [42]. As reviewed in Appendix A, these 
symmetries act on A^-point correlation functions in Fourier space as 

/ N \ 



N 

5k^M^ = i (2(A, - 3)% + kl^,^ - 2ka ■ 4%) A = . (2.13) 



a=l 



These constrain the form that the correlation functions take in momentum space. Here we assumed 
that the free evolution (2.11) dominates at late times. If this is not the case, as we discuss in Sec. 6, 
one cannot trade the time dependence of correlation functions for A's. 



3 Non-linearly realized conformal symmetry 

Returning to the action (2.7), we will be led to consider a particular choice of mass, corresponding 
to the quadratic action 

= Ml I d^xV^ ("2 + 2^'^') ' (^-l) 

which has = — 4i7^ (A = —1). In fact, this action is the quadratic action for the Goldstone 
of broken conformal symmetry. We briefly review how (3.1) can arise naturally on an effective 
de Sitter space through the spontaneous breaking of conformal symmetry in a flat space quantum 
field theory. We will consider two simple examples, and see that they both lead to the quadratic 
action (3.1) for fluctuations. 
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Negative quartic potential. As a first example, we consider (/>^ field theory on flat space 
with 'wrong-sign' coupling constant [19, 21, 29, 33]: 

S = I d'x l^-^idcpf + ^^'^ , (3.2) 

where A > 0. This action is classically invariant under conformal transformations,^ consist- 
ing of space-time translations P^, Lorentz transformations J^'^ , space-time dilation D, and 
special conformal transformations K^^. These 15 symmetries form the algebra SO (4, 2). 

Assuming a homogeneous field profile, the equation of motion for (j) reduces to </> — Xcp'^ = 0, 
with zero-energy solution 

Here H will soon be understood as playing the role of a fake Hubble constant. This back- 
ground, which describes the field starting out from the top of the potential in the asymptotic 
past and rolling down subsequently, spontaneously breaks 5 of the original conformal sym- 
metries, namely P^, J*^* and K^. The 10 unbroken symmetries, D, P*, J*-', and K^, form the 
SO {4, 1) de Sitter subalgebra. Expanding about the 1/t background as cj) = (p + we obtain 
the following quadratic action for fluctuations: 

d'x ^-lid^f + . (3.4) 

Through a field redefinition, (p = (j) + (f = 4>e'^ = nf-t) ^'^'^ introducing an effective de 
Sitter metric ^ 

3% = Jp^''^!^'^ ' (^•^) 



the quadratic action becomes 



which is precisely of the form (3.1). As advocated, this action arises through the spontaneous 
symmetry breaking 50(4, 2) — 50(4, 1). 

• Galilean Genesis. In its simplest guise, Galilean Genesis [20] is achieved with a (wrong- 
sign) kinetic term plus a cubic conformal galileon term: 

s = Jd^x (/2e™(an)2 + |Jnn(an)2 + £^idu)'^ , (3.7) 

where the scales /, A have dimensions of mass, and the scalar field H is dimensionless. This 
action is also invariant under the conformal group 50(4,2), but in this case dilations and 
special conformal transformations act non-linearly to start with. The equation of motion 
following from this action admits a background solution of the form [20] 

1 2A'^ 

This solution preserves the de Sitter subgroup of the conformal group. Perturbing about this 
solution H = rH-vr, and introducing the effective de Sitter metric (3.5), the quadratic action 
takes exactly the form (3.6). 

^Incidentally, this theory is also asymptotically free. 
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More generally, the action for vr follows solely from the pattern of symmetry breaking. In 
[33], the most general action for the goldstone of the breaking pattern 50(4,2) — t- S'0(4, 1) was 
constructed. To cubic order in the fields and second order in derivatives, it is given by (see Appendix 
B for details of the construction) 

S^ = Ml j d^xV=5 (-^(^^)^ + 2^^^^ - <dTTf + AH'^t:^ , (3.9) 

which is consistent with (3.1) at quadratic order. This action non-linearly realizes time translations 
, boosts J*^* and the temporal component of special conformal transformations , which act as 

djg^ir = -J + tdiTT - XidtTT ; (3.10) 

-,2 

SkoTT = - {2tx''^^. - X^dt) TT . 

Notice that the — naively unstable — growing mode solution vr ~ 1/f of (3.1) corresponds merely 
to a non-linearly realized time translation, and is therefore harmless. 



4 SO (4,2) SO (4,1) consistency relation 

The next question is: how do these non-linearly realized symmetries act on correlation functions? 
In this Section we show that the non-linear realization of conformal symmetry constrains the form 
that correlation functions take in the limit that one of the vr external legs is taken to be very soft. 

4.1 Derivation of the consistency relation 

In the case of inflation in the decoupling limit, the isometries of de Sitter are spontaneously broken 
by the inflaton's time-dependent background to the subgroup of rotations and translations. As a 
consequence of this spontaneous breaking, there are specific relations between correlation functions 
of different order. In particular, the (A^ + l)-point correlation functions in the squeezed limit 
are related to the variation of the A^-point correlation functions under the broken symmetries 
(dilations and special conformal transformations). These relations go by the name of consistency 
relations [3, 5-9]. They are the Ward identities resulting from the non-linearly realized symmetries 
in the broken phase of the theory [10-12]. 

Our aim is to show that similar relations hold in the case of the nonlinearly-realized 50(4, 2) 
symmetries. We again expect that the squeezed limit of an (A + l)-point correlation function is 
related to the action of the broken generators on the A-point function. In this case the broken 
generators are time translations P^, boosts J^^, and the time component of a special conformal 
transformation K^; correspondingly, the consistency relation will contain three pieces. 

In what follows we are going to use the background- wave arguments developed in [3, 6, 8]. In 
order to derive the consistency relation, the crucial insight is that, due to non-linear realization of 
<S'0(4, 2), the effect of a long background mode of vr can be obtained by a coordinate transformation 

(</.(xi) . . . 0(xjv))^L = (Hxi) ■ ■ • Hxn)) , (4.1) 

where (/> is a generic scalar field of the theory. As we will argue shortly, the change of coordinates is 
related to the value of the long mode as vrL = —5t/t. We can use the previous equality to express 
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the {N + l)-point correlation function in terms of the A^-point function 

(7rL(x)(/)(3;i) . . . 0(xAr)) = {lTi^{x)d{(j){xi) . . . (t){xN))-Ki) = (7rL(x)(5(0(xi) . . . <i){xN))) . (4.2) 

This relation among correlation functions of different order in real space is the starting point of our 
derivation. By computing explicitly the variation of the A^-point function on the right-hand side 
and then going to momentum space, we will obtain the consistency relation. 

Before diving into the details of the full calculation, as an illustration of the method we consider 
the simplest case where we have only broken time translation. The coordinate transformation 

t^t' = t + ao (4.3) 

will induce the change in the perturbation around the background 7r'(x') = Tr{x) — ao/t, thus 
generating the homogeneous mode 

VTL = -^^ . (4.4) 
On the one hand, the (A^ + l)-point function satisfies 

(tTl^I . . . (Pn) = {T^h5{<Pl ■ ■ ■ <AAf)7rL> • (4.5) 

On the other hand, the A^-point function in the presence of a long mode is the same as the A^-point 
function without the long background but computed in the transformed coordinates: 

6{(j)i . . . (I)n)tti^ = ^{(t>i ■ ■ ■ (l>N)t~^t+ao = "-o^i^i ■ ■ ■ ^n) ■ (4.6) 
Using the relation (4.4) between ttl and ao, and substituting the result in (4.5), we finally obtain 

{ni^cpi . . . (Pn) = -(7rLvrL)t^(0i . . . (/>Ar) . (4.7) 

This relation is written in real space but going to momentum space is trivial: a homogeneous 
long mode translates to the zero- momentum limit. Thus the (A^ -|- l)-point function has to be 
evaluated in the squeezed configuration 

■ ■ ■ hJ<i-.o = -i^.-^^yt^^ih, ■ ■ ■ ' (4-8) 

where primes indicate that the factor of (27r)^(5(^ ki) has been removed from correlation functions. 
The previous expression has the same form as Maldacena's consistency relation for standard infla- 
tion [5, 6]. This is not a surprising result. In the case of the decoupling limit of standard inflation, 
Maldacena's consistency relation can be seen as a consequence of broken dilation invariance. The 
squeezed limit of the (A^ -|- l)-point function is related to the scale transformation of the A^-point 
function. In our case, where time translations are spontaneously broken, we have a similar results 
— the squeezed limit of the (A^ -|- l)-point function with a soft vr is related to the variation in time 
of the A^-point function. 

We are now ready to derive the most general consistency relation by applying the above method 
to the full 50(4, 2) — )• 50(4, 1) symmetry breaking pattern. Since we have 5 broken symmetries 
but only one Goldstone, the consistency relation will contain several terms of different order in g, 
each coming from a different broken generator. Our starting point is the most general infinitesimal 
conformal transformation of space-time, which can be written as 

t^t' = t + ao + axi + 2{b^x^)t - bQ{-t'^ + x'^) = t + St , 
Xih^x'i = Xi + tti - Cit + JijXj + 2{bf^x^)xi - bi{-t'^ + x'^) = Xi + 6xi . (4.9) 
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Notice that in 5t and Jx,-, we can neg lect all 0{t) and 0{t^) terms because these will give subdom- 
inant contributions in the late-time limit. The remaining terms in 5xi consist of time- independent 
translation, rotation, dilation and 'id special conformal transformation, all of which are linearly 
realized symmetries. We conclude that for the purpose of deriving the consistency relation the only 
relevant part of the coordinate transformation is the time transformation: 

t^t' = t + aQ + axi - box"^ = t + 6t . (4.10) 

This induces a long mode ttl of the form 

T^hix) = -y = -^(«o + CiXi - box'^) . (4.11) 

Two important points are worth emphasizing about the shape of the long mode. First, unlike 
the case when only the time-translation is broken, the long mode is in general x-dependent. Three 
different broken generators correspond to three terms with different powers of Xi . At the level of the 
consistency relation, this will translate to different broken generators constraining different powers 
of q in the g — )• expansion of correlation functions. Second, the above ttl does not correspond to 
the most general shape of a mode expanded to second order 

7r(x) = ^(0) + diTr{0)x' + ^didj7r{0) x'x^ + ... (4.12) 

Specifically, our coordinate transformation induces only an isotropic profile for the long mode at 
quadratic order in coordinates (proportional to x^), while the traceless part of didjir^O) cannot be 
generated through the coordinate transformation. Consequently, the consistency relation can only 
inform us about quadratic terms that are averaged over all angles. We focus only on the trace part 
of the long mode 

^l(x) '= ^l(O) + di7TU0)x' + ldW{0)x^ , (4.13) 

o 

where here ^= means that we have averaged over angles for the terms which are second order in 
gradients. Comparing this profile with (4.11) we can read off the values of oq, Cj and bo: 

ao = -tvrL(O) ; a = -td^7Tj^{0) ; bo = Itd^UO) . (4.14) 

D 

We now have at our disposal all the necessary ingredients to derive the full consistency rela- 
tion. The first step consists of computing the variation of an A'^-point function under the coordi- 
nate transformation (4.10). Since the only relevant part of the full conformal transformation is a 
spatially-inhomogeneous shift in time, it is natural to first calculate unequal-time correlation func- 
tions and then take the equal-time limit at the end of the derivation. As another technical step, 
we must choose a point around which to perform our coordinate transformation. Since the fields 
in the correlation function are evaluated at xi, . . . , xat, a natural choice is their average location: 

N 



a=l 



l,E^a. (4.15) 



In this way the origin of coordinate transformation is translated in the region of space where the 
correlators are calculated. The choice of x+ is of course arbitrary, and will drop out of the final 
result. The upshot is that the coordinate transformation of interest is 

t^t' = t + 4 - c+(x+ - Xi) - 6+(x+ - xf . (4.16) 
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where the coefficients are of the form (4.14), but now evaluated at x+i clq = —t7:i^{x^), cf = 
—tdi7ri,{x+) and 6q" = t(9^7rL(x+)/6. 

In order to be consistent, we must average all quadratic terms — as a result, the 2bQX ■ x+ term 
will give zero contribution. The variation of an arbitrary A^-point function under this coordinate 
transformation is 

5{4){ti,Xi) . . .(j){tN,XN))^-^ =^ ^ (a^ - cf{xf - Xai) - &o (^+ + ^l))-^{Hh,Xl) . . . (t){tN,XN)) ■ 

a=l 

(4.17) 

Rewriting this equation in momentum space and using (4.14), we obtain 

where P = ki + . . . + is the total momentum and A is the Fourier space n-point function. 
Expanding the phase e*^+ ''+ in small momentum g+ and averaging all second order terms, we get 

(l - ilt4 - + iqfxai - \qlxi^ e«+-^^+ =^ 1 + iqfxai - \qlxl + 0{ql) . (4.19) 

The result is independent of x+, as anticipated. 

The final step in the derivation is to average over the long mode Tri^{t, x). After performing the 
average, integrating over g+ and replacing the coordinates in front of the exponent with derivatives 
with respect to the momenta, we obtain the {N + l)-point correlation function 

{n{t, xmti,xi) . . . <P{tN, x^)) =' / ^ ■ • • 1^ (03 (2vr)353(p)p^(^) 

X E (-*^'^(*^' • • • '^^)) - + l"^'^'^) e"^'^'^''-'" ■ (4.20) 

Next we integrate by parts to move the dk^^ derivatives from the exponent to the amplitude A. 
Since the delta function S^{P) depends on the momenta ki, . . . ,kN, this will introduce new terms 
involving derivatives of the delta function, but these are in fact precisely what is needed to build 
up the Taylor expansion of 6^{P + q). Indeed, averaging all second order terms, it is easy to prove 
the following relation^ 

(^1 + q^dk^^ + ^-q'dl^ [6HP)A') 6\P + gl (l + q^d,^^ + ^-q^dl^j A' . (4.21) 

Going to momentum space on the left-hand side of (4.20) and comparing the integrands on each 
side of the equation, we obtain 

i-<rmsh) . . . c/>,^^(t;v))Uo -tp^Q) E (i + ^^^^^^ + l^i'^l) wJ'f'kSti) ■ ■ ■ hjt^)y ■ 

(4.22) 

^Notice a change of the sign in the second term in the brackets coming from integrating by parts. 
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As expected, the (A^ + l)-point correlation function in the squeezed Hmit is related to the action 
of broken generators of the conformal algebra on the A-point function. Notice that on both sides 
the same delta function 6^{P + q) has been removed from correlation functions, which implies in 
particular that the N momenta on the right-hand side do not form a closed polygon. As we will 
see later with explicit checks, this prescription is essential for the consistency relation to hold. 
From (4.22), it is easy to deduce the equal-time version of the consistency relation: 

• • • ^kJ,-.o =' -PAq) (l + E + E 9l) 4t{h. ■ ■ ■ ■ (4-23) 

\ a a / 

This is the main result of the paper. In the remainder of this Section, we will verify the consistency 
relation in a few examples. Explicit calculations of the correlation functions will be given in 
Appendix C. 



4.2 Explicit checks of the consistency relation 

The simplest checks involve the squeezed limit of a three-point function. As a general rule, note 
that the three-point function in the squeezed limit does not have any linear correction in q. Indeed, 
using the symmetric expansion of ki and k2, 

ki = k-^; k2 = -k-^; ki + k2 + q = 0, (4.24) 

it is easy to see that (4.23) does not give a linear term in q. Consider the three-point function of a 
TT with two massless spectator field x given by (C.26). Its expansion in the soft limit, 

(-?(^)x.^(^)x,^,(t))U = 1^1^^ {^(k . q)- - k-e + 0(,3)) , (4.25) 

starts at order g^. Notice that the angular average of the 0{q'^) terms vanishes, in agreement with 
the consistency relation (4.23) — since the two-point function (x^^X^j)' is independent of time for 
massless fields, the right-hand side of (4.23) vanishes in this case. 

As a less trivial check, consider the three-point function of tt with massive (A = 1) fields 
given by (C.19). Its squeezed limit after averaging is 

/ V 9i/4 1 , 3(fc-gl^-fcV ^ avg 9H^ 1 

^^m,n,U^Q- 2MlMU''k\ 4A;4 )' 2M^M^q^k- ^ ' 

It is straightforward to verify that the consistency relation is satisfied at zeroth order in q using 
the explicit expressions for the power spectra given by (C.12) and (C.17). Similarly to the massless 
case, the C(g^) terms average to zero. On the other hand, from the right-hand side of (4.23) we 
conclude that the consistency relation seemingly predicts non- vanishing 0{q^) corrections. This is 
an illusion: remembering that the hard momenta do not form a closed polygon, ki and k2 are not 
equal in magnitude and opposite, and using (4.24) we obtain 

ki,2 ~ k i 2^2 8A:4 I • ^^-^'^ 
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Averaging, this gives 

1 avg 1 / ^ 1 o 1 avg 1 / 1 



12^ ^ k\ 24A:2 

Thus, we see that all 0{q'^) contributions to the right-hand side of (4.23) cancel, as they should. 

Similar relations can be used to check the consistency relation for the three-point function of 
the Goldstone field tt, given by (C.14). Its squeezed limit (after averaging) gives 

Meanwhile, expanding the right-hand side of (4.23) in the soft momentum q gives 

The consistency relation is satisfied. 

We can also check the consistency relation for the squeezed limit of the four-point function (C.22) 
with one Goldstone and three massive (A = 1) fields ^p. In this case we can express one of the 
momenta, say ^3, in terms of the others: /cs = —{ki + ^2) — 9- Averaging all terms proportional to 
q^, we obtain the following set of relations 

^3 p\ ) ' \ki + ^ h\ kl ) ^ 

where p = ki + k2- Using these equations, the four-point function can be expanded to 0{q'^) as 

avg SlTTH^Xt 1 ( h-q k2-q 2p-q 2q^\ 
^^<fPAn,n,)<i~^o 8M^M^ q^kik2P [ Skj Skj Sp^^dp^' ^ ' 

Meanwhile, from the right-hand side of the consistency relation we are instructed to act with various 
derivatives on the three-point function for massive fields (see (C.18)): 

Noting that the action of the Laplacian on monopole terms gives zero, d? = 0, and taking into 
account the power spectrum of the long mode vr, the right-hand side of the consistency relation 
becomes 

rhsof(4 23) - - ^^"^'^^ ' A _ i^j- _ _ 

^ ^ ^ 8M2M4 q^kik2k3 [ 3kj 3kl 3kl J 

Slvrif^At 1 / ki-q k2 ■ q 2p ■ q 2q 



avg 



8M2M4 q^kik2P r 3kl 3kj 3p2 9 



+ 7^h ■ (4-34) 



where in the last step we have used the expansion of k^. This verifies the consistency relation in 
this case. 
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Figure 2: Factorization of an (A^ + M)-point function via an exchange of a soft tt. 

5 Connection to observables: soft internal lines and 
anisotropy of the power spectrum 

As we discussed, the breaking of 50(4, 2) implies the existence of the Goldstone field vr, and con- 
sequently the consistency relation we derived constrains correlation functions with soft external 
Goldstone lines. Unfortunately, the cosmological perturbations we observe come from a spectator 
field and not from vr, so that it is not obvious how one can connect the previous results to obser- 
vations. There are, however, two situations in which 50(4,2) is observationally relevant. The first 
is when diagrams of the spectator field contain soft internal vr lines (for similar results in inflation 
see [3, 43, 44]). Internal soft vr lines are expected to give the dominant contribution when a sum 
of external momenta becomes small, and they will dominate in comparison with soft internal lines 
of the spectator field, because of the very red spectrum of the Goldstone. The second possibility 
stems from the fact that, even if vr is not directly measured, its value during the conformal phase 
is correlated with the modes of the spectator field and thus changes their statistics. In particu- 
lar, very long modes of vr induce an anisotropy in the spectator field power spectrum. These two 
observational features were studied in [37-39]. Here we want to stress that these properties are a 
direct consequence of the non-linear realization of 50(4, 2) and not specific to a given model. We 
will also find an additional important contribution to the four-point function from a loop of vr fields 
that has been overlooked in the literature. This contribution may be larger than the tree-level vr 
exchange and it is phenomenologically quite different. 

Let us start with soft internal vr lines. In the limit in which the sum of N external momenta 
becomes small, the amplitude of an (A^ + M)-point function factorizes in the following way (see 
Fig. 2) ^ 

^^fci • • • ^fcAZ+JV^'J-^O = ^^--f^fci • • • ^fcM^'J->o(^gX^^^^j • • • Xfc^,^^^)qr->0 • (5-1) 

The (A^ + 1) and (M + l)-point functions are severely constrained by the 50(4, 1) symmetry and 
their squeezed limit is further constrained by the non-linear realization of 50(4, 2). In this way, the 
amplitude for the (A^ + M)-point function with a soft internal line can be expressed in terms and 
M-point functions. The simplest case is the four-point function of massless spectator fields, which 
was studied in detail in [37-39]. Using the factorization (5.1) above and the squeezed limit (4.25) 
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for the three-point function (vrxx); in the hmit ki + k2 = we get 

(X£. . . . Xj;Uo = ^ (3(^. . <,f - l) (3(i3 . a? - l) . (5^2) 

It is important to stress that the shape of the four-point function in the soft internal hmit is 
completely specified by symmetries since the three-point function {nxx) is completely fixed by 
50(4, 1) up to an overall constant. Notice that the squeezed limit of the three-point function is 
constrained, as we discussed in the previous Section, by 50(4, 2) as well. In the massless case we 
cannot obtain terms scaling as q^ or q^ , and all terms scaling as q^ must vanish when averaged over 
the angles. This is indeed what we have in (5.2). 

The four-point function becomes very large in the q ^ limit, as it scales as 1/q. This is a 
consequence of the very red spectrum of tt and it can be contrasted, for example, with inflationary 
models with reduced speed of sound which are regular in the g — )• limit [45, 46] . We conclude 
that a four-point function which becomes large in the soft internal (collapsed) limit, with the precise 
shape (5.2), is a general prediction of the conformal mechanism. Notice, however, that the overall 
multiplicative constant in (5.2) cannot be fixed by symmetry arguments. 

The four-point function we studied is obtained by averaging over the long wavelength modes of vr. 
However, if we do not take the statistical average, we still have a realization-dependent effect: long 
modes of vr induce an anisotropy in the power spectrum of the short modes, as pointed out in [37- 
39]. Notice that this is possible even though tt does not contribute to the observed perturbations: 
its value during the conformal phase still affects the observable modes of the spectator field. This 
effect is also completely fixed, up to an overall constant, by the symmetries of the problem. 

The effect of a long vr mode on the observable 2-point function can be read from the three-point 
function {ttxx): given by (4.25), in the squeezed limit 

(""qXk.XkX-^o = I^^-^x^5Pt^^^^ ^ - 1) , (5-3) 

where 9 is the relative angle between k and q and = 9H'^/2M^ ; = H'^ /2M'^ are the 
dimensionless power spectra. We can write the variation of the power spectrum of x in the presence 
of a given realization of the vr field in a schematic way as 

^ix,X_,y = ^^^-^ = te-.^)' • ^^^3cos^^ - Dt^^vr,. (5.4) 

All modes vr^- which are outside the present Hubble radius will contribute to the anisotropy of the x 
power spectrum. The typical size of the effect is given by the square root of the variance calculated 
by summing over all super Hubble modes 

{tVTr.^-qY = A r tV^, ~ • (5.5) 



This gives 



(2vr)3 ' ^ 2^2 ^ . ^5^2 4^2 



= (XfcX^fe)' ( 1 + civ^^(3cos2 0- 1) ) , (5.6) 



where ci is a number of order unity, which depends on our position in the Universe [37-39]. 

Another source of anisotropy in the power spectrum arises by considering a four-point function 
(vrvrxx) [37-39]. This induces a variation of the 2-point function (xx) in the presence of two long 
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modes of tt. In this case the 5*0(4, 2) symmetry fixes both the shape and the normahzation of the 
effect. The variation of the 2-point function (xx) in the presence of two long background modes vri 
and TT2 corresponds to the composition of the associated 5*0(4, 2) transformations. A possible issue 
is that the broken generators Kq, Jqi and Pq do not commute, so that the overall transformation 
seems to depend on the ordering. Fortunately, in our case all the commutators of the broken 
generators of 5*0(4, 2) give unbroken generators. These do not change the 2-point function, so that 
we do not have to worry about non-commutativity in the case at hand. Since the 2-point function 
is time-independent, its variation at lowest order in gradients will come from a boost at second 
order. Without loss of generality, we consider a boost along the j;-direction. The transformation 
of coordinates is given by 

x = -t{x - Vxt) , y' = y, z' = z, t' = -f(t- v^^x) , (5.7) 

where 7 = (1 — w^)"^/^. Neglecting parts proportional to t, the induced background field is, up to 
second order in Vx, 

vr = -'^ = ^. (5.8) 
t t ^ ^ 

In momentum space, the parameter Vx is given by 

Vx = itqxT^q ■ (5.9) 

The transformation (5.7) implies that in momentum space kx component of the wavevector has to 
be multiplied by 7^^, while ky and k^ remain the same. Expanding fc^^ in the denominator of the 
power spectrum, we find that the effect on the 2-point function of x is: 

= (X,X.,)^^ = -{x,xy^tV4-os^ . (5.10) 
We can calculate the typical value of t'^q'^Tr'i- in a way similar to before: 

(*'?^a) = ^/ 9 dgt g ^P.log— , (5.11) 



(27r)3 ' ^ ^' 27r2 Jq ^ ^5^2 27r2' " A 

where A is an IR cutoff. The contribution to the anisotropy is given by: 

^(Xgx_,^)' = -(Xgx_g)^7^.1og§^cos20 . (5.12) 
Combining with (5.6), the total anisotropy of the power spectrum is given by 



l + ciVK-^(3cos2e-l) + C2Kcos2^1og-^ , (5.13) 
k Ay 

where C2 is another constant of order unity, which depends on the particular position in the Universe. 
The two sources of anisotropics are quite different. The first scales as 1/fc, and thus important only 
for long modes, while the second is scale invariant. Moreover, the first contribution averages to zero 
if summed over the possible orientations between long and short modes, while the second does not. 
Notice also that the first effect is dominated by vr modes which are slightly longer than the present 
Hubble radius, while the second gets contributions from all scales as shown by the logarithmic 
dependence. The logarithmic enhancement can overcome the suppression due to the fact that the 
second effect is of order tt^ and not tt. 
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Figure 3: The four-point function of x's with an exchange of one and two soft vr's. 



As we have seen, the second contribution to the power-spectrum anisotropy is related to the 
correlator (vrvrxx)- This suggests that we missed a potentially large contribution to the four-point 
function of x's in the soft internal limit, coming from a loop of soft vr particles (see Fig. 3). At 
first this looks worrisome as we expect a loop diagram to be small compared to a tree-level one. 
However, the situation is similar to the one we discussed for the anisotropy. When only one soft 
vr is exchanged, the interaction with the x's arises at order as we discussed above. When two 
soft vr's are exchanged, on the other hand, each of them carries a single soft momentum, as the 
interaction arises from the non-linear realization of boosts. Therefore, in going from tree-level vr 
exchange to a one-loop diagram the number of g's at the vertices remains the same, and we have 
the extra loop factor 

(2vr)3 g5 ~ g2 • ^^-^^^ 

If q is small enough compared with the external momenta, the loop diagram will dominate over 
the tree level exchange. Notice that this does not signal a breakdown of perturbation theory: it is 
straightforward to check that the exchange of extra vr's is not further enhanced by powers of 
but only suppressed by powers of P^. 

The loop diagram is straightforward to evaluate starting from (5.10) 



Akf J (2vr) 



(5.15) 



where q = ki + k2 and + ^2 = Q- In writing this expression we have assumed that both internal 
legs are soft so that their coupling is fixed by the non-linear realization of 50(4, 2). Indeed we will 
see that the loop integral is dominated by having qi and q2 both of order q. If we disregard the 
angular dependence and average over the direction of the short modes, we get 



4 kf kl 



d=^gi 
(2vr)3 



[Ql -92) 

Qi Q2 



(5.16) 



The loop integral is dominated by long modes and it is IR divergent, similarly to what happened 
for the anisotropy of the power spectrum. We get 



(5.17) 



As promised this result contains, when compared with the tree- level calculation (5.2), a factor of 
Vnk'^/q'^ which may be large for sufficiently small q. 



17 



Notice that the momentum dependence of this result (after performing the angular average) is 
exactly the one of a tnl non-Gaussianity, with 



p2 




(5.18) 



This four-point function will therefore contribute both to a stochastic scale-dependent bias [40] and 
to the power spectrum of ^-distortion [41]. It would be interesting to understand whether the angle 
dependence, which is different from a standard tnl shape, affects these observables. 



We now address two important subtleties which arise in our analysis. The first (Sec. 6.1) is that, 
due to the negative mass of the Goldstone field vr, it is prima facie unclear when it should be 
possible to consistently assign 3d conformal weights for the transformations (2.13). We will see 
that this is only possible when the relevant interaction turns off sufficiently quickly, so that the 
in-in integral is dominated by early-time contributions. 

The second subtlety (Sec. 6.2) arises when checking the consistency relation: we must expand 
the correlation function slightly off-shelF to check the relations, but it is not clear that there is a 
unique way to do this. We show that not only is the procedure unambiguous, but it is possible to 
work on-shell at all times by rewriting the consistency relation at each order in g as a differential 
identity. These identities correspond precisely to the Ward identities derived from the broken 
symmetry generators [47]. 

6.1 Assigning conformal weights 

We first discuss a puzzle that arises when considering correlation functions involving the Goldstone 
field IT. It should have conformal weight A = — 1. Indeed, the correlation function {ttxx) with x 
massless satisfies (2.13) with the proper A's. On the other hand this is not the case for (vrvrTr) or 
{iTip(p) with = 1. These correlators are 50(4, 1) invariant but one is not allowed to replace the 
time dependence with the free-field conformal weights. 

The resolution of this apparent paradox is as follows: the possibility of assigning conformal 
weights to each field as a function of its mass depends on whether the late-time behavior of cor- 
relation functions is fixed by the free-field evolution. In other words, interactions must turn off 
sufficiently fast in time. Since one expects in general that the late-time evolution is purely classical, 
it is straightforward to ascertain by inspection whether conformal weights can be assigned. 

To make this a bit more concrete, consider the three-point function 



where -ffj^/ is the cubic, interaction-picture Hamiltonian. The dependence on t comes both from 
the integration limit and from the (j) fields. This second contribution follows the free-field evolution 
and allows a consistent assignment of conformal weights. Thus, we want to study under which 
conditions the additional t dependence from the integral can be neglected, i.e., whether the integral 

^Hcrc, our usage of the phrase "ofT-shell" is slightly nonstandard. What we mean is that the momenta in 
the correlation functions do not form a closed polygon. 



6 Subtleties 




■oo(l — ie) 




(6.1) 
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is dominated by early times. Notice that the ie prescription is different in the two complex- 
conjugated terms, but we can neglect this difference as we are only interested in the late-time 
behavior of the integral. Hence we can write 

{Mt)Mt)Mt)) ~ -^ fdt' {[Mt)Mt)Mt),Hl^ht')]) 

= -i fdt' {Ut)ct>2{m,it),H[^{t')])-i fdt' {Mt)[Mt),H^!ht')]Mt)) (6.2) 

-i fdt'{[Mt),Hl^l{t')]Mt)Mt)). 

In each of the three terms, one leg is evolved non-linearly, and we want to show that this describes 
the classical non-linear evolution. Let us focus on the first term for concreteness. The commutator 
with the Hamiltonian will contain three pieces, depending on which of the three legs ^3 is paired 
with. For each of these contributions, we can integrate by parts to get rid of the derivatives acting 
on the leg we singled out and write the Hamiltonian as 

hI^ = jd'xaU{t')0^'\t'), (6.3) 

where 0^'^\t') is a quadratic operator which will be paired with the legs 4)i and </>2. In the late-time 
limit, we can replace 4>i and <j)2 within 0^'^\t') by their growing mode solutions, which we take to 
be real. Thus 0^'^\t') becomes a number which acts like a source for the mode cj)^. We have the 
mode expansion 

0(t') = 0el(t')«^ + 4>ll{i')^ > (6.4) 

with (pel denoting the properly normalized wavefunction. Expanding the commutator, we obtain 
-i fdt'{[Ut),l d'xa^<pO^^\t')]) = -i fdt'a'[Mt')€iit)-Mt)€i{t')]O^^Ht'). (6.5) 

Inside the brackets we recognize the Green's function, so that this expression gives exactly the 
classical non-linear evolution of the field (j)^. In the late time limit this Green's function has two 
contributions, of the form 

The first contribution depends on t as the growing mode solution of the free theory.^ Therefore, 
this contribution will give the time dependence which follows from the assignment of conformal 
weights provided that the integral in (6.5), namely 

fdt'j-/^^0^'\t'), (6.7) 

converges. This, as we discussed, means physically that the classical interaction is dominated by 
early and not late times. Notice that for = and 0^^\t') = const, we have a logarithmic 
divergence. A positive m? or an interaction containing derivatives makes the integral converge. 
What about the second contribution (~ t^+t'^-) in the Green's function? This is always suppressed 
compared with the one we discussed as \t'\ < \t\, and they become comparable at late times t' ^ t. 
Again, if the above integral converges, this late time limit is irrelevant, and the second contribution 
to the Green's function is negligible. 

what follows we assume that rr? jH.'^ < 9/4 so that A± are real. 
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6.2 Off-shell ambiguity 

As we have already pointed out, the momenta on the right-hand side of the consistency relation 
do not form a closed polygon. In other words, the correlation function has to be evaluated slightly 
"off-shell" . However, when calculating the A^-point function, there are many different ways one can 
go off-shell.^ To see this, let us write a generic A-point function in the following form 

((/>fc^ • • • </>fc^) = (27r)35(^i + . . . + %)^' . (6.8) 

Adding a term that is zero on-shell to the amplitude M does not change the result. For an arbitrary 
function F, we can write 

(<A,-^ . . . 4>i;i = i27Tf S\P) (^A' + P- F{h, ...,kN)) , (6.9) 

where, as before, P = hi. However, once we go off-shell, it is far from clear that the right-hand 
side of the consistency relation gives the same result independently on the choice of the function 
F. This has to be the case for the consistency relation to make sense, because the left-hand side is 
unambiguous. Below we will prove that there is in fact no ambiguity in the consistency relation, 
and we will show how it can be rewritten in a form where all the correlation functions are calculated 
"on-sheh" . 

The simplest check is to explicitly compute the action of the differential operator in the con- 
sistency relation on an expression that is zero on-shell and that can be generically written in the 
form P ■ F{ki). A straightforward calculation gives 

\ a a / 

= iP + q).F+^^^Yl ^'^.-^^ + ^9^^' H ^^^^^^ + ^^'P^ E • (6.10) 

a a a 

Using the fact that P -|- g = 0, we see that the first term on the right-hand side is zero, while the 
second and third cancel upon averaging. The last term is of order 0{q^) and can be neglected. 
Indeed, as we expected, adding to the amplitude of an A-point function a term that vanishes on- 
shell does not change the result on the right-hand side of the consistency relation (up to corrections 
of order 0{q^)). This proves that there is no ambiguity, and that different choices of going slightly 
off-shell give the same contribution in the end. 

We are now going to give a slightly more general proof of the same statement, which will make 
more explicit the constraints that different broken generators impose on the relation among correla- 
tion functions. The aim is to start from the consistency relation written in form (4.22) and rewrite 
it as a set of different relations that constrain different powers of q. If we understand expression 
(4.22) as a Taylor expansion of the (A -|- l)-point function, then this procedure corresponds to 
isolating the coefficients in this expansion. Since we have constraints on momenta from the delta 
function on both sides, and we are computing one of the correlation functions slightly off-shell, this 
procedure is far from trivial. As a first step, let us rewrite the consistency relation (suppressing 
for simplicity the corresponding time dependence of the modes that is the same as before) in the 
following form 




(6.11) 



'We thank Alberto Nicolis for bringing this ambiguity to our attention and for discussions on the matter. 
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At zeroth order in q this relation trivially becomes 

1 



(6.12) 



where both sides are calculated on-shell. As expected, this expression is the same as (4.8). The 
right-hand side can be understood as a zeroth order in a Taylor expansion and, obviously, in this 
case there is no ambiguity in the consistency relation. To find the next coefficient in the Taylor 
expansion of the left-hand side, we must differentiate with respect to q and set q = 0. 



da 



1 



N 



q=0 



a=l 



(6.13) 



q=0 



There are two things to be noticed. Firstly, using the on-shell condition ki + . . . + g = we 
can write dq- = — 9^^^.. Secondly, since we set q to zero, the prime on the correlation function on 
the right-hand side now means that this correlator is calculated on-shell. The previous expression 
becomes 



da 



1 



7V-1 



g=0 



a=l 



(6.14) 



The difference of the derivatives in the brackets can be rewritten in a form of a total derivative 

d^ATi „ d 



kNi 



dk 



-dkMj 



dka 



(6.15) 



Since the correlator is on-shell and any dependence on k^ can be removed, we can trivially include 
into the sum a total derivative with respect to /c^r. The final expression that we obtain is 



da 



1 



N 



q=0 



a=l 



(6.16) 



The last relation represents the consistency relation for broken boosts, and it corresponds to the 
terms in (4.23) that are first order in q. The most important point to stress is that, in this form, 
the correlation functions on both sides are calculated on-shell, and there is therefore no ambiguity 
in the final result. Finally, we can follow the same procedure to isolate the terms of order g^. The 
resulting expression is 



1 



P.{q) 



N 



q=0 



a=l 



(6.17) 



Equations (6.12), (6.16) and (6.17) are equivalent to the consistency relation (4.23). As we have 
already pointed out, in this form all correlation functions are computed on-shell. This is another 
proof that in the consistency relation (4.23) there is no ambiguity in off-shell prescription on the 
right side of the equation. These derivative identities are in fact the Ward identities corresponding 
to the non-linearly realized time translations, boosts and special conformal transformation, respec- 
tively. A formal derivation of these relations using the machinery of [11] will appear in future 
work [47]. 

To perform explicit checks using the derivative form of the consistency relation, we work on- 
shell on both sides at all times. Therefore, on the left-hand side, we express one of the momenta. 
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say kN, as a sum of the other momenta. We then take the squeezed Hmit, obtaining a function of 
q, the small momentum. In order to check the various relations, wc can then take derivatives of 
this left hand side with respect to q and then set q = 0. On the right-hand side, we must also work 
on-shell. This means that we also write the /cjv momentum in terms of the other N — 1 momenta 
(not including q). 

Schematically, the procedure is as follows: consider checking the consistency relation 



a=l 



d d 



ih. ■■■hJ' ■ 



(6.18) 



Wc rewrite the left hand side so that it is a function of N different momenta, that is we take 
kjsi = —J2^a — Q- We then take the squeezed limit q ^ and differentiate with respect to q. On 
the right hand side, we write k]\r = —J2^a- This means that we actually only have to take N — 1 
derivatives on the right hand side. 

For illustrative purposes, we provide an explicit check of the consistency relation in differential 
form. Consider the soft limit of the three-point function involving only tt fields, (tt'^). The three 
and two-point correlation functions are given by 



81H^ (g^ + kl + fci) 
4M4 q^klklt^ 



9H^ 1 



(6.19) 



We take the limit g — )• to obtain the squeezed limit of the three-point function 

9iJ2 2 5(q-ki) ^[n<i-k?-q^ki) 



2 , 5(g-fci) 

1.5 h.7 



2k\ 



(6.20) 



From this, we immediately read off: 

1 



vr 



_d_ 
dq^ 



1 



PAq) 
1 

pM 



q=0 



g=0 



m'^ 1 

H^kft^ " 
45i72 k{ 

90^2 1 



-t 



1 d d 

2 dfcfdt 
d 



1 d2 

t- 



(6.21) 



M2 klt^ 2 dkjdV ^1 



where the last step in each equation follows from differentiating {tt^^tt^J' . Thus the derivative 
form of the consistency relation checks out at each order. Note that when checking the identities 
in Ward form, it is not necessary to take an angular average in the 0{q'^) terms. The Laplacian, 
9q, appearing on the left hand side automatically picks out the trace part of the order q'^ terms of 
the soft mode. 



7 Conclusions 

Symmetry considerations allow us to make powerful, model-independent statements. In this paper, 

wc have considered how non-linearly realized symmetries in the conformal mechanism enforce var- 
ious relations among the correlation functions of the theory. In doing so, wc have uncovered some 
model-independent and robust predictions of the mechanism. It is worthwhile summarizing them: 
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• Absence of detectable gravitational waves. 

• Model-dependent local non-Gaussianity from the conversion mechanism. 

• Anisotropy of the power spectrum, see (5.13), [37-39]. 

• 4-point function in the soft internal limit due to tree-level vr exchange, (5.2), [37-39]. This is 
relevant on large scales. 

• 4-point function in the soft internal limit due to one-loop vr exchange, (5.15). This dominates 
for sufficiently small internal momentum, and it shows up as stochastic bias and in the power 
spectrum of fi distortion. 

It will be interesting to study which of the last three signatures is the most constraining, as they 
all depend on the same parameters. See [48] . 

It will also be interesting to further elucidate the properties of correlation functions when more 
than one external leg is taken to be soft. These higher-soft limits can serve to probe the underlying 
broken symmetry algebra, and could provide another set of consistency relations, not just for the 
conformal mechanism, but also in the case of inflation. 

Another interesting direction to pursue is to derive the consistency relations for the DBI version 
of the conformal mechanism studied in [27, 28]. In this scenario, the symmetry breaking pattern 
is the same, but the non-linearly realized conformal symmetries act in a different way. It would be 
interesting to understand if the consistency relations are the same in this case, and whether the 
generic predictions of the conformal mechanism, such as statistical anisotropy, go through in this 
case. 
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A Conformal transformations on correlation functions 

In this Appendix we derive the action in Fourier space of the linearly-realized dilation and spatial 
special conformal transformations on correlation functions. 

A.l Dilation 

We will work in an arbitrary number of dimensions, d. The dilation 
in position space as 

6d(I> = (A - x^dA)<p . 
We note that the field (p can be written in Fourier space using 

we may therefore write 



operator acts linearly on fields 

(A.l) 

(A.2) 
(A.3) 
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Now, we can integrate by parts to obtain two terms 

5d^ = J d'^ke'^''' (^A - d - k ■ dk^ 4>k 
Prom this, we deduce the Fourier space transformation rule 

SD(t>k = - (^{d- A) + k- 4) 4>k 



(A.4) 



(A.5) 



Now, we want to obtain the action of dilation on a correlation function. A correlation function has 
two parts, the amplitude and the delta function, schematically it is of the form 



6dA = Sd{s^P)A') , 



(A.6) 



where the prime indicates removal of the delta function and P is the sum of the momenta P = J2k. 
We may then write 



Sd {s^{P)A') =-Y.{(d- ^a) + ka ■ 4) [5\P)A!' 

a=l 

N 

= -A'P ■ dp5\P) - {{d - A„) + ka ■ 4) A' . 



a=l 



The term outside the sum may be integrated by parts to obtain a factor of d. The term where the 
derivative dp hits A' vanishes because it multiplies P5^{P) = 0. We then have 



r ^ 

5d {s^P)A') = 6'{P) d-J2{id- Aa) + K ■ 4.) A' 



a=l 



(A.7) 



Prom this, we deduce that the dilation operator acts on the amplitude without the delta function 
as 



SdA' 



N 

-d{N -l)+J2(Aa-ka- 4„) 



a=l 



A' 



(A.8) 



A. 2 Special conformal transformations 

Special conformal transformations act in real space as 

6ka = {2Ax^ - 2xAX^dB + x^dA)4> ■ 



(A.9) 



Following the same steps as above, we may write this in Fourier space acting on the primed correlator 
as 



N 

SkaA' = (2(A„ - d)d,A + k^Bi^ - 2ka ■ dkA^) A' 



(A.IO) 



a=l 
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B Construction of actions for spontaneously broken 
conformal symmetry 

In this section we review the construction of actions for the spontaneous breaking of the conformal 
group down to its de Sitter subgroup [33]. Equivalently, this construction gives actions for fields 
on de Sitter space that non-linearly realize conformal symmetry. The basic idea is fairly simple, 
in order to linearly realize the isometries of de Sitter, we consider a de Sitter metric g'^. A scalar 
field action constructed using this metric will be invariant under the de Sitter group. In order to 
non-linearly realize conformal symmetry, we add the conformal mode to this metric and consider 

9^.u = e^^g''^ . (B.l) 



B.l Action for the Goldstone 

In order to construct the action for the field vr — the Goldstone mode of broken conformal symmetry 
— we construct curvature invariants from the metric g^y. Since the metric is conformally flat, all 
of the information is contained in the Ricci tensor 

R^, = 3H^gf, - 2V^V,7r - gf^OTT + Id^ird^n - 2gf,{dT^f . (B.2) 

Then, in order to construct actions for vr, we merely form any diffeomorphism scalar using the 
conformal metric, its covariant derivative and the Ricci tensor. The lowest order lagrangian is 
given by the invariant volume and the kinetic term is given by the Ricci scalar 

~ y=5 (i? - m') = ^/=5ds (^e2-(9^)2 + K'-Ut. - H'e'- + ^e'^^ , (B.3) 

Where we have added the cosmological term so that there will be no tadpole about vr = 0. Ex- 
panding this lagrangian out to cubic order and integrating by parts reproduces (3.9) in the text. In 
[33], it was shown that this construction is equivalent to the coset construction of Callan, Coleman, 
Wess, Zumino and Volkov [30-32]. This construction misses a four-derivative Wess-Zumino term 
(see [49] for details). 



B.2 Coupling matter fields 

It is also straightforward to couple matter flelds to the Goldstone field vr; one thing we can do is 
use the covariant derivative of the conformal metric (B.l) to construct actions for matter fields 
(f), contracting indices with g^,^. Additionally, we are allowed to multiply curvature invariants 
constructed from R^i, by arbitrary polynomial functions of cf).^^ The two-derivative lagrangian for 
a scalar field is then of the schematic form 

£0 ~ ^ (-^(9<A)' - 1/(0) + /(</')i?) = (-^e'"(9</))' - e'^VicP) + e''^f{^)R^ . (B.4) 

Since our aim is to merely check the consistency relations in a variety of examples, we consider the 



m 



2 

case where V{(/)) = -^(j)'^ + Xc/)^ and /(</>) = 0. The lagrangian then takes the form 

(1 TfT^ \ 

-^e'^d<pf-^e'-^'-Xe'-cP^\ . (B.5) 

^°Note that since the WZ term shifts by a total derivative, we are not allowed to multiply it by a function. 
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Expanding about </) = vr = to quartic order yields the action 

S^ = MlJ d'^x./^ (^-lid^? - ^</'' - 2m2vr02 _ ^(dc^f _ Ac/>3 - AXttc/)'^ . (B.6) 

B.3 Transformation of tt 

It is also quite straightforward to work out the way the non-linearly realized conformal symmetries 
act on TT. Under an infinitesimal diffeomorphism, the metric changes by the Lie derivative 

Sgfiu = g^iu = -gpv^^ii'' - gt^p^ui'' ■ (B.7) 

We assume that the background metric g'^ remains fixed (this restricts us to isometrics of de Sitter 
plus conformal transformations), so we have 

25^g^, = -gpuVf^^P - g^^pV^^" , (B.8) 

tracing over both sides gives Jvr = — jVp^^. This is the divergence of a vector, so we may write 



Stt = — . d, 



p 



So we have the transformation rule for vr 



-ge) = -edpTT-\vfip. (B.9) 



Jvr = -edpTi - hjfe . (B.IO) 

From this transformation rule, it is clear that that vr will transform linearly under isometries of the 
dS metric (V^^^'' = 0) and will transform in a nonlinear fashion under broken transformations. To 
make this explicit, we must make a choice of de Sitter slicing. Choosing the planar inflationary 
slicing: 

we find 



5n = -(Pdp7T-^df,e + le ■ (B.12) 



From this, it is straightforward to derive the symmetries (3.11). 



C Correlation functions 

Here we collect some results for correlation functions involving spectator fields coupled to the 
Goldstone field vr. 



C.l Mode functions for massive fields 

In this Appendix, we derive the expression for the mode functions of a massive scalar field on de 
Sitter space in terms of Hankel functions. This expressions are needed to compute the correlation 
functions we need to check the consistency relations. Consider the general quadratic action for a 
massive scalar 
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Where is an arbitrary mass. The equation of motion fohowing from this action is 



2 ■ ■Til 

We define the canonically-normalized variable 

Ma , 



Hi-t) 



(C.2) 



(C.3) 



whose mode functions satisfy 



k'- 2 



i/2 t2 



Vk = . 



(C.4) 



Defining x = —kt and 
Bessel's equation 



I — , after changing variables to fi^ = v^/ ^/x this can be cast as 



(C.5) 



which is well-known to be solved by (we choose Hankel functions as our basis) 

fk(x) = ci{k)Hi^\x)+C2{k)Hl^\x) . 



(C.6) 



We fix the coefficients by demanding that in the far past {—kt — t- oo), the mode functions of 
the canonically normalized variable, Vk, have their Minkowski space form. This is the so-called 
adiabatic vacuum (Bunch-Davies) choice. That is, we demand 



-ikt 



(C.7) 



Then, using the asymptotic expansion for the Hankel functions as — fct — t- oo 

^ 1 



^(i)(_A;t)^_ef(i-)y| 



-ikt 



yf-ki 



Akt 



vr yj-kt 



This implies that we need to take ci{k) = — e 2 (2 and C2{k) = in (C.6). This leads to 

the expression for the (j)k mode functions 



-e 2 I2 



where Hi^\-kt) is a Hankel function of the first kind. Note that for m > the solution is a 
Hankel function of imaginary order. 



'9 "1^ 

4:~lP 

9H^ 



(C.8) 
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C.2 In- in integrals 

In order to compute correlation functions, we employ the Schwinger-Keldysh or in-in formalism 
(see [5, 50] for an exposition). In this formalism, rather than computing in-out S-matrix elements, 
we compute correlation functions sandwiched between the same vaccum. The correlation function 
for an operator, 0{t) is given by [5, 50] 

{0{t)) = (O|re'^*o'^*'^-'(*')o(t)re"'^*o'^*'^'"*^*'^|0) . (C.9) 

Here Hi^t is the interaction Hamiltonian, T denotes time-ordering while T denotes anti-time- 
ordering and to is an early time. Generally we will only work to leading order (tree-level) where 
the correlation function is given by 

{0{t)) = -i f dt' (0 I [0(t),/7int(t')] I o> . (C.IO) 



C.3 Correlation functions of tt 



Here we compute the two and three-point correlators for the Goldstone field vr. We consider the 
action (3.9). The quadratic equations of motion lead to the following mode function for the field tt 



vrfc(t) 



V 4 5/2^ 



-kt) 



1 + ikt 



-3H 

/W{-t)M^ 

From this the two-point function can straightforwardly be computed: 

9H^ 



2+2 



kH 



-ikt 



P-n{k) = {tTj:tt_ 



2MlkH'^ 



kh^ kH"" 



Note that this field has an extremely red spectrum, peaking strongly as /c — )■ 0. 

From the action (3.9), we can also compute the three-point function, (vr^) 
Hamiltonian, -fTintj at this order is minus the lagrangian 



int 



d^x^i- 



Mi / d'^x 



1 



TT{dTr) 



TT 



Applying the formula (C.IO), we obtain (at late times) 

81H^ {kf + kl + fcf ) 



4M4 kfk^k^t^ 
where is a cutoff introduced to regulate the divergence as t — t- 0. 



(C.ll) 



(C.12) 



The interaction 



(C.13) 



(C.14) 



C.4 Massive spectator field, A = 1 

The simplest case of a spectator field coupled to vr is a massive field with = = 2H^, 
corresponding to 3d conformal weight A = 1. We take the action (B.6) with this choice of mass: 

S^ = Ml I d^xV^ [-\{dv? - - 4i?V(/52 - ^{d^f - Xif^ - AXinpA , (C.15) 
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The mode functions for the field are given by 
which leads to the two-point function for the spectator 



We can also compute various higher-point correlation functions involving this spectator. The 
simplest is the three-point function involving only the tree-level correlation function is given by 

Additionally, we can compute the (tt^^p) three-point function for these fields. There are two con- 
tributions to the correlation function, one from each of the -mpip vertex and the 7r(5(/?)^ vertex; the 
final result is given by 

This correlation function is invariant under (Ad) dilations and under 5j^i with Aq = { — 1, 1, 1}, 
agreeing with our general arguments for when conformal weights may be consistently defined, in 
spite of the fact that this correlation function does not scale in the naive way with time. 

Finally, we compute a four-point function, involving three (/? fields and one Goldstone; this 
computation is slightly more involved. There are two contributions to this four-point function, one 
coming from a contact diagram involving the nip^ vertex and one coming from an exchange diagram 
at second order in the vertices involving a single vr and two (/j's. The interaction Hamiltonian is 
given by^^ 

= Ml J (^T^') . (C.20) 

The correlation function is then a sum of three terms 

/t 
-oo 

+ r dt' f dt"(0|F(,^,)(t')vr,-v^,-,¥',^,V',^3(t)^£^(i")|0> (C.21) 

J — oo J —oo 

-2Re(f dt' f dt''(0k,-v.,^^c^,^^^,^3(t)F(f)(t0^S(O|0)) . 

\J—oo J —oo J 

When the dust settles, the four-point function is given by 



^^Note that at this order, we must be careful in deriving the interaction Hamiltonian, in this case it is still 
minus the interaction lagrangian, but in general this will not be true at quartic order. 
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c.5 Massless spectator field, A = 

We now consider a massless spectator field, corresponding to (B.6) with with = A = 0. The 
cubic action for this field is given by 

S^ = M^J d^x^ - <dxf^ . (C.23) 

The mode functions for x are the well-known result for massless fields 

Using this, the two point function for a massless field is the standard result 

Pxik) - (X£X_g)' = ^^(1 + kh') . (C.25) 

Additionally, we can compute the three-point function {ttxx) using the standard techniques, sum- 
marized above. At late times, we obtain 



q^{kf + k'l) - {kf + kl) + 3{kfk'l + kjk^) 



9H^ 1 
~ 8M2M2 q^kfk^ 

This correlation function is invariant under 6j^i with = { — 1,0,0}. Additionally, it has the 
leading scaling behavior with respect to time that is expected. 
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